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ABSTRACT 

The emission line ratios [Olll] A5007/H/3 and [Nil] A6584/Ha have been adopted as an empirical 
way to distinguish between the fundamentally different mechanisms of ionization in emission-line 
galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal 
structure of the emitting Hll regions, and these calculations depend on the assumptions one makes 
about the relative importance of radiation pressure and stellar winds. In this paper we construct 
a grid of quasi-static Hll region models to explore how choices about these parameters alter Hll 
regions' emission line ratios. We find that, when radiation pressure is included in our models, Hll 
regions reach a saturation point beyond which further increases in the luminosity of the driving stars 
does not produce any further increase in effective ionization parameter, and thus does not yield any 
further alteration in an Hll region's line ratio. We also show that, if stellar winds are assumed to 
be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is 
inconsistent to simultaneously assume that H II regions are wind-blown bubbles and that they have 
high ionization parameters; some popular Hll region models suffer from this inconsistency. Our 
work in this paper provides a foundation for a companion paper in which we embed the model grids 
we compute here within a population synthesis code that enables us to compute the integrated line 
emission from galactic populations of Hll regions. 

Subject headings: galaxies: high-redshift — galaxies: ISM — HII regions — ISM: bubbles — ISM: 
lines and bands 



1. INTRODUCTION 

The line ratios [Olll] A5007/H/3 and [Nil] A658 4/Ha 
first proposed fo r use in galaxy classification by |Bald-| 
win et al. (1981) (hereafter BPT), are commonly used 
to diagnose the origins of emission lines from galaxies, 
and in particular to discriminate between galaxies whose 
emission is powered by star formation-driven H II regions 
and from those powered by active galactic nuclei (AGNs). 
These emission line pairs are particularly useful because 

(1) they are bright and thus relatively easy to measure, 

(2) blending between the lines can be corrected with rea- 
sonable accuracy, so long as the spectra are taken with 
sufficient resolution, and (3) the wavelengths in each line 
pair are quite similar, so the line ratio is relatively inse n- 
sitive to dust-reddening ( Veilleux fc Osterbrock|[l987 l. 

The power of these line ratios as diagnostics comes 
from their sensitivity to the spectral shape of the radi- 
ation field driving the ionization, which can be under- 
stood from a simple physical picture. To first order, 
the intensities of the Ha and H/3 lines simply measure 
the total photoionization rate, and thus normalize out 
the ionizing luminosity. On the other hand, the [Olll] 
and [Nil] intensities are sensitive not only to the total 
ionizing luminosity, but also to the shape of the ioniz- 
ing spectrum and to the ionization parameter U, which 
measures the ratio of photons to baryons in the ionized 
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gas. When the ionizing flux arises from hot stars, the 
ionizing spectrum is dominated by low-energy photons 
that have short mean-free paths through neutral gas. 
Thus the H II region consists of a fully ionized zone with 
a sharp boundary. Within this region, as IA increases, 
more of the ionized gas volume becomes filled with high 
ionization-potential species such as ++ , and less with 
low-ionization potential species such as N + . As a result, 
Hll regions ionized by hot stars tend to fall along a se- 
quence that runs from high [O III] /H/3, low [Nil] /Ha to 
low [O III] /H/3, high [Nll]/Ha. On the other hand, if 
the ionizing spectrum follows a power-law, as expected 
for AGN, then a significant amount of the ionization is 
produced by X-ray photons capable of ionizing higher 
ionization potential species like ++ . Moreover, these 
photons have large mean-free paths, giving rise to a large 
zone of partial ionization rather than a smaller region of 
full ionization as in the stellar case. In this configuration, 
[O III] /H/3 and [Nil] /Ha both increase with Li, and ei- 
ther one or the other tends to be larger than in the stellar 
case, leading to a sequence that runs from intermediate 
to high [O III] /H/3 and [Nil] /Ha and is well-separated 
from the locus occupied by H II regions dominated by 
stellar sources. 

This simple picture is roughly consistent with local 
observations: star-forming galaxies in the Sloan Digi- 
tal Sky Survey (SDSS; median redshift z = 0.1) obey a 
tight correlation between the [O III] /H/3 and [Nil] /Ha 



ratios in the BPT di agram (Brinchmann et ajJ 2004 



Tremonti et al. |2004 ). However, higher redshitt star 



forming galaxies are offset from this sequence to higher 
values of [OlIl]/H/3, without joining the locus of points 
occupied by AGN in the SDSS sample (Shapley et al. 



20051 JErb et al.| |2006l [20101 |Brinchmann et aP 2008 ; Liu 
ct al. 2008pSs our work is motivated by a desire to bet- 



wind gas is of ten smaller rather than larger than radia- 



ter understand the physical information encoded in the 
BPT diagram, we pause to consider how this shift in the 
BPT locus might arise. 

The difference in line luminosity ratios could be intrin- 
sic to the galaxies' Hll regions. As we have said, the 
ionization parameter is a major controlling factor which 
positions regions along the star-forming locus. Line emis- 
sion from these regions will be also affected by the metal- 
licity and dust content of the interstellar gas; by the 
density of that gas (through the critical densities of the 
lines); by the ionizing spectra of the stars (which reflect 
stellar masses, metallicities and rotation rates). In ad- 
dition, stellar winds can alter the bou ndary conditions 
for ion ized zones, a point discussed by |Yeh fc Matzner| 
(2012 hereafter YM12) and to which we return below. 

Alternatively, the shift in the BPT diagram could arise 
from outside the H II regions if their light is mixed with 
line emission f rom shoc ks or an unresolved AGN (e.g., 
Liu et al |2008| . Indeed, [Wright et al.| ( [2010| use integral 



held spectroscopy to demonstrate that a weak AGN is re- 
sponsible for the sh ift in a single galaxy at z = 1.6, and 



Trump et al. ( 2011 ) stack HST grism data of many galax- 
ies to show that this phenomenon is reasonably common. 
Taken together, these studies raise the possibility that 
Hll regions at z ~ 2 lie along the same BPT locus as 
those nearby, and the shift is an optical illusion caused 
by active nuclei. 

However, the distribution of high-redshift galaxies in 
the BPT diagram is also shifted in the direction of high 
U. Because radiation pressure rises, relative to gas pres- 
sure, in proportion to U, this implies that the radia- 
tion force typically is more important in high-redshift 
galaxies. This radiation-force-dominated condition is 
also more prevalent among starburst galaxies in the lo- 
cal Universe, as YM12 argue on the basis of mid-infrared 
line emission. This possibility has also received signifi- 
cant support from recent resolved observations of H II re- 
gions, which provide direct evidence that radiatio n pres- 
sure is signi ficant for the most luminous examp les (iLopcz 
et~aTJ|2l)lTl but also se e |Pellegrini et^[20TTl and |ljilich 



fc Tenorio- tagle|[2013 ) . 

The detailed role of radiation pressure in altering the 
line ratios of starlight- ionized Hll regions has received 
relatively little attention, although the phenomenon has 



(Binette et al.11997 


Dopita et al. 2002 ) . Early models ig- 


nored radiation pressure entirely (e.g.|Dopita et al. 200U 



diation pressure (e.g. Dopita et al.|2005 


2006a b 


Groves 


et al. 


2008 


Levesque et al. 2010p, it is 


either explicitly 



or implicit 



y assumed that the geometries and internal 



structures of Hll regions are dominated by stellar wind 
bubbles rather than radiation pressure. As we discuss 
below, the assumption that stellar wind pressure exceeds 
radiation pressure is often physically inconsistent with 
the range of ionization parameters being probed. More- 
over, resolved observations of the brightest nearby Hll 
regions indicate the hot gas produced by shocked stel- 
lar winds for the most part does not remain confined 
within Hll r egions, and instead leaks out into the low- 
den sity IS M ( |Townsley et al.|2003| [Harper-Clark fc Mur- 
|ray||2009) . As a result, the pressure of shocked stellar 



tion pressure ( Lopez et al.||2011 



YM12 



In this paper we explore how radiation pressure influ- 
ences the line emission of H II regions. To do so, we com- 
pute a sequence of hydrostatic H II region models under 
a variety of physical assumptions about the relative im- 
portance of radiation pressure and stellar wind pressure 
(Q, and we explore how varying the physical assump- 
tions alters the loci occupied by the model H II regions in 
the BPT diagram (©. We then compare our models to 
those published by other authors (Q and draw conclu- 
sions (Sj5j). In a companion paper (Verdolini et al. 2012, 
hereafter Paper II) , we use the grid of H II region models 
presented in this paper to construct a population syn- 
thesis model capable of predicting the line ratios of star- 
forming galaxies containing many different H II regions. 
We use these models to compare to observations of star- 
forming galaxies. Although this study cannot replace a 
full investigation of the factors affecting the H II region 
locus within the BPT diagram, it is the first to explore 
the roles of radiation and wind pressure. 

2. PHOTOIONIZATION MODELS 

2.1. Input Parameters and Calculations 

To study the influence of radiation pressure on H II re- 
gions, we construct a grid of static, single Hll regions, 
with a wide range of sizes, ionizing luminosities, and wind 
strengths. Our procedure is as follows. We first use 



the stellar popula tion synthesis code Starburst99 (Lei- 
therer et al. 1999) to generate spectra from coeval star 
clusters. We assume that all of the clusters are massive 
enough to fully sample the stellar initial mass function, 
which we take to have exponents -1.3 and -2.3 between 
stellar mass boundaries 0.1, 0.5, and 120 M Q . We em- 
ployed the Geneva standard evolutionary tracks (|Char- 
bonnel et al.|19 96 ; Schacrcr et al.|199 3a b ; Schallcr 'et al. 
1992) with solar mctallicity, and Lejeune-ScTrinutz stella r 
atmospheres (Lejeune et al. 1997 1998; Sch mutz||1998[ ), 
which incorporate plane-parallel atmospheres and stars 
with strong winds. We record the Starburst99 output 
spectra for cluster ages of to 11 Myr at 0.5 Myr inter- 
vals. 

We then use t he photoionizat i on co de Cloudy 08.00, 
last described by |Ferland et al.] ( |1998[ ), to compute the 
structure of static, spherical HII regions driven by point 
sources whose spectra are taken from the Starburst99 
calculations. In addition to the spectrum of the driving 
source, Cloudy requires a number of other input param- 
eters. The first of these is the total luminosity of the 
ionizing source, for which we run a series of models with 



L 



10 



33 



10 erg s in 1 dex steps. The second 



is the number density of hydrogen nuclei at the inner- 
most zone of the H II region, which we set to values from 
UH.in = 10 _1 — 10 5 in steps of 1 dex. The third is the dis- 
tance of the innermost zone from the point source, which 



we vary from R m = 10 2 f c h 



— 10 2 f c h in steps of 0.2 dex. 



4 Note that |PelIegrini et ah] |201l[ | assume a smaller filling 
factor for the X-r ay emitting gas, and so assign it a much higher 
pressure than |Lopez et aL] p011). For the same luminosity, small, 
higher-pressure bubbles have a greater dynamical effect on their 
immediate surroundings. But these bubbles are less important 
for the entire region than large, lower-pressure ones. This is a 
consequence of the virial theorem, which ties dynamics to the net 
energy budget. 
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Table 1 

Comparison of Model Parameters 





Dopita et al. 2000 


Model Our Model 


Levesque et al.|2010 


Model 


Ionizing Spectra 


Starburst99 


Starburst99 


Starburst99 




Initial Mass Function 


Salpcter, M up = 120 M Default Starburst99 IMF a , M up = 120 Mg 


Salpeter, M up = 100 Mn 


Stellar Tracks 


Geneva standard Geneva standard 


Geneva high mass 


-loss 


Stellar Atmospheres 


Lejeune-Schmutz Lejeune-Schmutz 


Pauldrach/Hillier 


Photoionization Code 


MAPPINGS III Cloudy v08.00 


MAPPINGS III 


Radiation Pressure 


Not included Optional 


Included 




Geometry 


Plane-Parallel Spherical 


Plane-Parallel 


H II Region State 


Uniform gas pressure Hydrostatic 


Isobaric b 




Metallicity 


Solar 


Solar 


Solar 




Gas Phase Abundances (logarithmic, relative to H) 


H 













He 


-1.01 


-1.01 


-1.01 




C 


-3.74 


-3.74 


-3.70 




N 


-4.17 


-4.17 


-4.22 







-3.29 


-3.29 


-3.29 




Ne 


-3.91 


-3.91 


-3.91 




Mg 


-5.12 


- 5.12 


-5.12 




Si 


-5.45 


-5.45 


-5.45 




S 


-4.79 


-4.79 


-5.19 




Ar 


-5.44 


-5.44 


-5.44 




Ca 


-8.16 


-8.16 


-8.16 




Fe 


-6.33 


-6.33 


-6.33 





a IMF exponents: 1.3 and 2.3 at mass boundaries 0.1, 0.5, and 120Mq. 

k L10 included radiation pressure in their isobaric calculations, thus the hydrostatic condition reduces to 
pressure in plane symmetry (YM12). 



state of uniform total 



Here the characteristic radius 



aL 2 



Teh 



12k (2. 2k B Tc) 2 S 



(1) 



is the radius of a uniform-density, dust-free Stromgren 
sphere for which the gas pressure is equal to the unat- 
tenuated radia tion pressure at its edge (Krumholz & 
Matzncr 2009; YM12); in this equation, a is the recom- 
bination rate coefficient, T is the gas temperature, and 
L and S are the bolometric luminosity and the output 
of ionizing photons per unit time, respectively, from the 
point source. In all the calculations presented here we 
adopt the same fiducial parameters as YM12: T = 8000 
K, a = 3.0 x 10~ 13 cm 3 s _1 . Finally, we adopt Cloudy's 
default ISM dust grain abundance and size distributions 
at solar metallicity, but in order to allow more meaningful 
comparison between our results and those of D00, we ad- 
just the gas-phase element abundances in our calculation 
to match theirs. The se choices me an that the dust dis- 
criminant parameter ( Draine|2011 ) takes the same value, 
7 = 7.4, as in YM12: so, dust opacity is significant within 
radiation pressure-dominated ionized zones. We summa- 
rize all the parameters we use in our calculations in Ta- 
ble |T] The Tab le also describes the p arameter choices 
usedin D00 and Levesque et al. (2010 hereafter L10). 

For each set of input parameters, we use Cloudy to 
calculate the structure of the resulting H II region, halt- 
ing at the point where the gas temperature drops to 100 
K in order to ensure that the ionization front is fully 
enclosed. We perform each calculation twice, once with 
radiation pressure turned off, and once with radiation 
pressure turned on and allowed to exceed gas pressure 
(in contrast to Cloudy's default setting, which does not 
allow radiation pressure to be greater than gas pressure.) 

2.2. Model Outputs and Physical Parameters 



The output of our calculations is two four- 
dimensional grids of models defined by the parameters 
(t, "H,in, L, i?i n ), where t is the age of the stellar popu- 
lation used to generate the ionizing spectrum. One grid 
contains models with radiation pressure turned on, which 
we refer to as the RP models, and the other contains 
models with radiation pressure disabled, which we re- 
fer to as gas pressure, or GP, models. For each run in 
both model grids, we compute several optical emission 
line luminosities integrated over the ionized gas, includ- 
ing the lines used to construct the BPT diagram: Ha, 
H/3, [Om]A5007, and [Nll]A6584. 

In order to understand the physical meaning of the 
results, it is helpful to characterize each model by two 
dimensionless numbers that can be computed from the 
Cloudy output. Following YM12, we define the radiation 
pressure parameter 

(2) 

where Rip is the radius of the ionization front (IF). A 
value of ^ < 1 indicates that the entire IF falls within 
f c h, and thus that radiation pressure is more important 
than gas pressure in determining its structure. Again 
following YM12, we define a separate stellar wind pa- 
rameter 

P V 

q = 111 111 fV) 

PwV lF -P in V^ {6) 

where Pif and P; n are the gas pressures at the edge of the 
ionization front and the innermost zone, respectively, and 
Vif = (4/3)7ri? 3 F and V in = (4/3)7i\R 3 n are the volumes 
contained within the IF and the inner edge of the Hll 
region, respectively. The inner edge of the photoionized 
region is the outer edge of the bubble of hot gas inflated 
by the stars' winds. Thus fl reflects the contribution 
of a pressurized wind bubble to the total energy budget 
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Figure 1. An example of how our models fill the parameter space 
of ^ and Q. In the Figure, each plot symbol shows the values of 
and Q computed for a particular calculation in our model grid 
at t = 0, nH,in = 10 cm -3 , with radiation pressure on. Colors 
indicate lines of constant L, running from 10 33 - 10 46 erg s 
as indicated in the legend. The sequence of points along a given 
model corresponds to varying _Ri n from 10~ 2 f c h — 10 2 f c h, with 
increasing with _Ri n . Note that if is a function of density as well as 
ionizing source luminosity. Radiation pressure can be significant 
in a high density region with relatively lower ionizing luminosity. 
Note that our full HII region model grids are incorporated into 
the dynamical models in Paper II. 

of the Hll region. In a region strongly pressurized by 
stellar winds !1 > 1, while in a region with negligible 
wind pressure Jl <C 1. 

For each Cloudy model, we compute the quantities "J 
and fi, and thus we may think of our models as describing 
a parameter space (t, nu t - ln , O), as illustrated by Figure 
[T] This parameter space describes H II regions for which 
both radiation and wind pressure run from strong to neg- 
ligible. To study the effects of winds, we reduce this four- 
dimensional parameter space to a three-dimensional one 
by selecting two representative values of f2: we designate 
models with log £1 = 2 as strong wind (SW) models, and 
those with log ft = —1.5 as weak wind (WW) models. 
Since our models never produce log £7 = —1.5 or 2 ex- 
actly, we construct these models by interpolation. At 
each age t, density nn.in, and luminosity L, we find the 
two models whose values of £1 bracket our target one, 
and we compute line luminosities at the target value of 
£7 by interpolating between the two bracketing models. 

Through this procedure, we obtain a set of four re- 
duced model grids, which we refer to as RPWW (radia- 
tion pressure turned on, log £1 = —1.5), RPSW (radiation 
pressure turned on, log $1 = 2), GPWW (radiation pres- 
sure turned off, log J! = —1.5), and GPSW (radiation 
pressure turned off, log£l = 2). 

It is important to bear in mind that our RP models are 
physically self-consistent, whereas the GP models are de- 
liberately not. Thus, RPWW models make a transition 
from classical spherical Stromgren spheres to ra diation- 
confin ed shells, along the sequence described by |Draine| 
(2011 1, as \& decreases through unity. RPSW models are 
always thin shells: both radiation and wind pressure play 
a role in confining them, but, as we explain below, wind 
pressure always dominates. GPSW models are also thin 
shells, but due to the neglect of radiation pressure they 
sample a range of ionization parameters inaccessible to 
real regions. Finally, GPWW models are always filled 



Stromgren spheres, even when the radiation force should 
confine them. They can also sample unphysically high 
values of U. Our GP models have strictly uniform gas 
pressure, as they include no other forces. 

Each of these model grids gives the line luminosities of 
H II regions as a function of the three remaining parame- 
ters, (t, riH,in> or equivalently (t,nn,L). We summa- 
rize the properties of the models in Table[2| We will make 
use of the four full model grids in Paper II, but for the 
remainder of this paper we concentrate on the particular 
case t = 0, nH,in — 10 cm -3 , in order to understand how 
the choice of input physics alters the structure of Hll 
regions. We choose these parameters in particular be- 
cause they match the ones used by a number of previous 
authors, thus facilitating easy comparison. 

3. RESULTS 
3.1. One-Zone Models 

For a given spectral shape, each ionized parcel with 
uniform density and temperature can be characterized 
by only two parameters: the density hh and the ioniza- 
tion parameter IA = n-y^/nu, where ?T"y^2 IS the number 
density of ionizing photons. Therefore there is a unique 
mapping between IA and initial densities on the BPT dia- 
gram. The one-zone models also represent simple analogs 
of H II regions, for one can decompose an H II region into 
zones in which U, Uhj and ionizing spectrum are nearly 
constant. Thus the one-zone models represent thin, uni- 
formly ionized regions which are very much like the ion- 
ized layer of a wind bubble. As such, they resemble best 
the SW models to be discussed in §3.2.2. 

We run an additional set of Cloudy "one-zone" cal- 
culations in which we only compute the properties of 
line emission from the first, innermost zone. In this 
zone we can specify the value of U by choosing the 
the density iih and the bolometric luminosity L (and 
thus the ionizing photon luminosity S). We run mod- 



els with nn = 10 



10 cm 3 in 1 dex steps, and 



U = 10~ 4 — lO 6 in 0.2 dex steps, all using an input 
spectrum corresponding to our t — Starburst99 model, 
and using the same abundances and other parameters as 
the rest of our models. 

In Figure [2] we show the constant U and constant nn 
contours marked with black solids lines and blue dashed 
lines, respectively, on the BPT diagram computed with 
the one-zone models. The ionization parameter U and 
ionizing luminosity S increase from lower right to upper 
left, and increasing the density shifts models up and to 
the right, until the density exceeds ~ 10 4 — 10 5 cm~ 3 . 
Beyond this point, the models shift down and to the left, 
because the density exceeded the critical densities of the 
[ N II ] and [ O III ] emission lines, whi ch are 6.6 x 10 4 cm~ 3 
and 6.8 x 10 5 cm -3 , respectively ( |Osterbrock fc Ferland 
2006 1. In Paper II, we will return to the discussion ot line 



ratios, U, and critical densities on the BPT diagram, and 
further discuss most extreme H II regions exceed i ng th e 



upper limit of line ratios set by Kewley et al. 
which is based on the mapping between line ra" 



(2001) 



ios and 



U but ignored the effect of densities. 

We note that the BPT locations of macroscopic Hll 
regions will differ from those of individual gas parcels, 
because of spatial variations in the physical quantities. 
In all cases U drops and the ionizing spectrum changes 
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Table 2 

Model Properties 



Model 


Acronym 


Radiation Pressure? 


log ft 


Radiation pressure, weak wind 


RPWW 


yes 


-1.5 


Radiation pressure, strong wind 


RPSW 


yes 


2 


Gas pressure, weak wind 


GPWW 


no 


-1.5 


Gas pressure, strong wind 


GPSW 


no 


2 




-0.5 

log 10 [NH]/Ha 

Figure 2. Models in the BPT diagram. Black lines show one- 
zone models with constant U, while blue dashed lines show one- 
zone models of constant njj; both are calculated for an ionizing 
spectrum corresponding to a zero-age stellar population and Solar 
metallicity. We also show models RPWW, RPSW, GPWW, and 
GPSW (orange and green lines, as indicated by the legend) with 
"H,in = 10 cm -3 , calculated with the same ionizing spectrum and 
metallicity. 

as one approaches the ionization front, because of selec- 
tive absorption by neutral H atoms, and in some cases by 
dust grains. When radiation pressure is strong (and is 
included) and winds are weak, the gas density increases 
significantly across the layer. We therefore anticipate 
that full Hll region models should differ from the one- 
zone calculation, even though the innermost zones are 
accounted by it. Moreover the macroscopic physical pa- 
rameters, the assumed geometry, and the inclusion or ne- 
glect of radiation pressure should affect the BPT loci. We 
explore these dependencies in the subsequent sections. 

3.2. Full Models 

We now turn to our four full (radially resolved) Hll 
region models, RPWW, GPWW, RPSW, and GPSW. 
In Figure [2j we overlay these models with t = and 
"■H.in = 10 cm~ 3 on the one-zone calculations. Other 
choices of density give qualitatively similar results, as 
long as the density is well below the critical densities of 
the [ O III ] and [ N II ] lines. As with the one-zone models, 
the full models form a sequence of values defined by 4" or 
IA, which we control by varying L: high-^, low-W, low-L 
models are found at the bottom right and low- 'I', high-W, 
high-L ones at the top left of each sequence. 

3.2.1. Weak Wind Models 

When stellar wind pressure is negligible (log SI = 
— 1.5), Hll regions at the low L end of the sequence are 
very similar to each other. This is because L determines 
the balance between radiation pressure and gas pressure; 
a high luminosity produces a large f c h (Equation [lj and 
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10' 
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RPWW, log £2 = -1.5 

GPWW, log £2 = -1.3 






1 

1 
1 



-RPSW, log £2 = 2.0 
GPSW, logQ = 1.8 



10' 
10 



10 

R 4 „ (cm) 



Figure 3. Electron density versus radius for sample HII regions. 
Top panel: RPWW (blue solid line) and GPWW (black dashed 
line) regions. Bottom panel: RPSW (blue solid line) and GPSW 
(black dashed line) regions. The age of the ionizing star cluster in 
these regions is Myr and the density at the inner boundary is 
10 cm" . The luminosity in all models is 10 43 erg s — 1 . We select 
the value of R ln from our grid that gives log Q closest to —1.5 and 
2; exact values of Q for the four cases shown are as indicated in 
the legend. See Section |2.2| for details. Again we note that high 
luminosity is required here to reach radiation pressure-dominated 
state because the density is low. 

thus a small value of \& (Equation[2 1 . Thus when L is low 
radiation pressure forces are negligible, and the results do 
not change much depending on whether we include them 
or not. The density within both Hll regions is roughly 
constant at % = «h in = 10 cm" 3 . 

At the high L, on the other hand, RPWW and GPWW 
differ substantially. In model GPWW, as L increases, we 
find that [OlIl]/H/3 increases and [Nil] /Ha decreases 
without limit. In contrast, in model RPWW these line 
ratios saturate at a finite value. If one were to infer 
ionization parameters from these line ratios based on 
one-zone models, one would say that U saturates at a 
finite value in model RPWW, while in model GPWW 
it can increase without limit as L does. We can under- 
stand the difference in behavior by examining the den- 
sity structures of RPWW and GPWW regions, of which 
we show an example in Figure [3| At high L, RPWW 
model H II regions are strongly dominated by radiation 
pressure. Under force balance, radiation pressure con- 
fines ionized gas into a much thinner layer and leads to 
a steep increase in density towards the IF. Much of the 
line emission comes from this dense layer, within which 
hi is much lower than it is closer to the central source. 
YM12 discuss this effect in detail. 

This effect does not operate in the GPWW models, 
where we have artificially disabled radiation pressure. As 
a result, these H II regions remain at nearly constant den- 
sity regardless of the source luminosity. This allows U to 
increase without limit, and in turn allows the [OlIl]/H/3 



G 



and [Nil] /Ha line ratios to continue changing even at 
large L. 

Finally, it is interesting to note that, despite the uni- 
form pressure in the GPWW models, the actual values of 
the line ratios are still significantly offset from the corre- 
sponding one-zone models of the same density, njj = 10 
cm~ 3 . At small L the shape of the sequence is similar 
but the models are displaced to slightly higher [Nil] /Ha 
and [OlIl]/H/3, while at large L the deviation is larger 
and the shape of the sequence is different as well. This 
difference occurs because, even though the pressure is 
uniform in the GPWW models, other quantities are not. 
In particular, the spectrum of the ionizing radiation field 
varies with radius, due to selective absorption of lower- 
energy photons by neutral H atoms and of higher-energy 
photons by dust within the ionized layer. 

3.2.2. Strong Wind Models 

In the RPSW and GPSW models (log 12 = 2), strong- 
stellar wind pressure produces large "voids" of diffuse, 
high temperature stellar wind gas at the centers of the 
model Hll regions. As a result, the ionized gas is con- 
fined to a thin shell between the wind bubble and the 
IF. 

The location of the RPSW model in the BPT dia- 
gram is strikingly far from the locations of other mod- 
els. Like the RPWW models, the RPSW models sat- 
urate at finite values of [Nll]/Ha and [OlIl]/H/3, re- 
gardless of how high the luminosity becomes. However, 
unlike in case RPWW, the saturation values are ex- 
traordinarily far down the sequence of one-zone models: 
[Nil] /Ha > 1CT - 5 and [OlIl]/H/3 < 10°, correspond- 
ing to a one-zone value of U < 10~ 33 . We can under- 
stand this effect by considering the relative importance 
of radiation and wind pressure in controlling the internal 
structures of Hll regions. A value of log 12 = 2 requires 
that PwVyp / PiuVi-a = 1.01. Physically, this amounts to 
saying that the energy of the wind bubble constitutes 
99% of the internal energy of the entire H II region. We 
note that V\f is strictly greater than l^n- Similarly, Pip 
is strictly greater than P n , since the radiation force nec- 
essarily falls to zero at the IF, and thus pressure balance 
requires that gas pressure at the IF exceed that at the 
edge of the wind bubble. Thus models with log 12 = 2 
necessarily have both Vxf ~ V- ln and Pif « Pi n . This cor- 
responds to the Hll region being a thin shell of nearly 
constant gas pressure. Figure [3] shows an example of this 
uniform density. 

The RPSW configuration clearly cannot have radiation 
pressure as a significant force. If the radiation pressure 
force were significant, then we could not have Pif ~ -Pin, 
since the pressure at the IF would be pure gas pressure, 
and this would have to balance the gas plus radiation 
pressure at the inner edge. The conclusion of this anal- 
ysis is that it is not possible to construct a physically 
consistent model in which wind pressure and radiation 
pressure are both strong in the dimensionlcss sense. In- 
deed, our model grids reflect this fact in that there are no 
models with radiation pressure turned on that are simul- 
taneously dominated by wind (12 3> 1) and dominated by 
radiation pressure (>& <C 1). This physical effect mani- 
fests in the BPT diagram as a saturation in the range of 
line ratios that the RPSW models are able to reach. As 
discussed above, the location of an Hll region driven by 



a stellar source in the BPT diagram is effectively con- 
trolled by U, the photon to baryon ratio. However, U 
may also be thought of as a measure of the importance 
of radiation pressure, since increasing the photon number 
density relative to the baryon density also increases the 
radiation pressure relative to gas pressure. The fact that 
RPSW models cannot reach small values of ^ also means 
that they cannot reach large values of and thus can- 
not reach the line ratios associated with large U. YM12 
used this point to derive upper limits on the wind energy 
budget within individual H II regions and entire galaxies. 

In contrast, radiation pressure is completely neglected 
in GPSW models. As there are no other forces to bal- 
ance gas pressure gradients, these models have Pif = 
P n independent of the luminosity, and thus a value of 
log 12 = 2 simply implies that the shell is very thin: 
Pif = 1.003P; n . In these models one can achieve ar- 
bitrarily high IA by raising the luminosity and increas- 
ing Pin to keep up with Pif- Figure [3] shows an exam- 
ple of such a model. Thus the GPSW models are not 
restricted in the range of U they can represent. How- 
ever, the comparison with the RPSW models shows that 
GPSW models at high U are unphysical, because radia- 
tion force would have compressed the gas and limited IA 
in a real region. 

The comparison between RPSW, GPSW, and one- 
zone models shown in Figure [2] also reflects these effects. 
Both RPSW and GPSW models are wind-dominated and 
therefore have nearly uniform densities, and thus can be 
reasonably-well represented by one-zone models, leaving 
aside the issues of radiation field hardeni ng an d non- 
uniform temperature discussed in Section 3.2.1| Thus 



both RPSW and GPSW models follow the one-zone se- 
quence reasonably closely. They differ only in the range 
of IA values within that sequence that they are able to 
sample. 

4. COMPARISON TO PREVIOUS WORK 

It is interesting to revisit earlier published models for 
Hll region line ratios in the context of our exploration 
of how these line ratios respond to changes in the in- 
cluded physics. We have chosen our model parameters 
in Starburst99 and Cloudy to be as close as possible to 
those used by D00 in order to facilitate this comparison. 
In Figure [4j we compare our four model results with the 
results computed by D00 and L10, at the same age of ion- 
izing star cluster (0 Myr) and same density nH.m = 10 
cm~ 3 . 

The D00 model, which did not include radiation pres- 
sure and adopted plane-parallel ionized gas slabs at a 
fixed density, is essentially a wind-dominated model. 
This is because a plane-parallel slab can be thought of as 
a thin shell of material at roughly fixed distance from the 
ionizing source, and the only way to create a thin shell 
of constant gas pressure is to confine it with hot gas. 
Therefore our closest model to the D00 model is GPSW, 
and indeed we find that our GPSW results agree with the 
D00 model fairly well. Differences in line ratios between 
D00 and GPSW are around 0.1 to 0.2 dex. Our model 
sequence extends somewhat further, but this is simply a 
result of our having used a slightly larger range of input 
luminosities. The agreement between the models con- 
firms that our Cloudy calculations, with input parame- 
ters set very close to the D00 settings, can reasonably 
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Figure 4. A comparison of model results on the BPT diagram. 
All results shown are for njj i n = 10 cm -3 and a spectrum corre- 
sponding to a zero-age stellar population. The D00 model (marked 
as KewleyOl in the legend) is shown in dark blue, the L10 model 
result is in light blue, and our model results are shown in orange 
and green lines. 

well reproduce the earlier results. However, we note the 
comparison shows that the D00 models are not physically 
realistic at high luminosity, because one cannot neglect 
radiation pressure in very bright Hll regions. Radiation 
pressure limits the physical range of U, particularly for 
wind-confined slabs, and models without radiation pres- 
sure such as those of D00 do not properly capture this 
effect. 

The models from L10 are similar to those of D00 in that 
they are based on plane-parallel ionized gas slabs, but the 
L10 models include radiation pressure. Therefore when 
U is low, the regions must be confined by wind pressure 
(like our RPSW). On the other hand, when U is high, 
the L10 models should be confined by radiation pressure 
(like our RPWW). However, Figure [4] shows that overall 
LlO's models closely track our GPSW curve (maximum 
separation < 0.1 dex). In light of our results, we can see 
that the L10 models, while not physically inconsistent, 
do represent a rather odd cut through parameter space. 
There are two structural parameters describing H II re- 
gions, and the L10 models sample a one-dimensional path 
through this two-dimensional space. Along this path the 
ratio of wind pressure to radiation pressure varies sys- 
tematically from large values (CI ^> 1, ^ ^> 1) at low 
ionization parameter to small values (!1 « 1, >P « 1) 
at high ionization parameter. There is no obvious phys- 
ical reason such a systematic variation in wind to radi- 
ation pressure strength should occur, particularly since 
the ratio of stellar wind moment um flux to luminosity is 
roughly the same for all O stars ( Repolust et al.p"004| ). 



5. CONCLUSIONS 

We have computed a grid of quasi-static Hll region 
models using Starburst 99 and Cloudy that covers a large 
range of density, luminosity, and stellar population age. 
In order to understand how radiation pressure and stellar 
winds alter Hll regions' internal structures and observ- 
able line emission, we run two sets of models, one with 
radiation pressure enabled and one with it disabled, and 
we vary the radius at which the inner, wind-dominated 
bubble ends and the photoionized region begins. In the 
manner, we construct four sets of model Hll regions: 



(1) ones with radiation pressure and weak stellar winds 
(RPWW), (2) wind bubbles that also include radiation 
pressure (RPSW), (3) Stromgren spheres where radiation 
pressure is ignored and winds are weak (GPWW), and 
(4) wind-dominated bubbles where radiation pressure is 
disabled (GPSW). We then explore how each set of Hll 
regions populates the BPT diagram. 

Our models reveal a number of interesting effects. 
All models form a sequence that runs from the lower 
right corner of the BPT diagram (high [Nil] /Ha, low 
[OlIll/HjS) to the upper left corner (low [Nil] /Ha, high 
[OlIl]/H/3), with the position of an Hll region along 
the sequence dictated by its luminosity, or equivalently 
its effective ionization parameter U. However, the range 
of U explored by the models is limited when radiation 
pressure is included. Because strong radiation pressure, 
which would produce high U, also causes gas to pile into 
a dense shell, the characteristic value of U within the 
shell is limited at a finite value. (See YM12 for more de- 
tail.) As a consequence, models which neglect radiation 
pressure can reflect an arbitrarily high value of U, which 
real regions cannot. 

The interaction of winds with radiation pressure fur- 
ther enhances this effect. We show that a stellar wind- 
dominated region cannot also have strong radiation pres- 
sure while remaining in hydrostatic balance, and as a 
result the range of 11 is severely limited. This means 
that wind-dominated Hll regions can never occupy the 
upper-left portion of the BPT diagram, and, conversely, 
those Hll regions that are observed to lie in this re- 
gion must either have negligible wind pressure, be far 
from pressure balance, or be kept in pressure balance by 
forces other than gas and radiation pressure (e.g. strong 
magnetic pressure; YM12). The most realistic option, 
and the one favored by direct observations of nearby 
Hll r egions (Harper-Clark & Murray 2009 Lopez et al 



2011 1 as well as mid-intrarared line ratios ( YM12) is the 
first one: wind pressure is not dynamically significant, 
at least for bright Hll regions. Further, the fact that 
the high-redshift galaxy population has characteristically 
high ionization parameters implies that radiation pres- 
sure is significant within these galaxies' ionized zones, in 
an ionization-weighted sense. 
We have compared our results to the earlier models 



of Dopita et al.| ( |2000[ D00) and |Levesque et al. ( |2010 
L10). In these models the Hll region is assumed to be 
a wind-dominated thin ionized shell, which corresponds 
to our GPSW model. We find that this model agrees 
well with the results of D00 and L10. However, we show 
that these models are inconsistent at the high luminosity 
end. The D00 models neglect radiation pressure for Hll 
regions where it is non- negligible. The L10 models in- 
clude radiation pressure, but we show that the assumed 
plane-parallel slab geometry is physically realistic only 
if the strength of the ratio of stellar wind pressure to 
radiation pressure varies systematically with H II region 
properties in a physically unexpected manner. 

While these calculations provide insight into how the 
physics driving Hll regions' structures translates into 
observable properties such as line ratios, a full model for 
where galaxies fall in the BPT diagram requires atten- 
tion to Hll regions' dynamical expansion as well as their 
internal structure. This problem is the subject of Paper 
II. 
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